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Proceedings contribution for "Dark Matter/Dark Energy 2004" 
ABSTRACT 

High-resolution, numerical simulations of 17 cluster-sized dark-matter haloes in eight different cosmologies 
with and without dynamical dark energy confirm the picture that core halo densities are imprinted early during 
their formation by the mean cosmological density. Quite independent of cosmology, halo concentrations have a 
log-normal distribution with a scatter of ~ 0.2 about the mean. We propose a simple scaling relation for halo 
concentrations in dark-energy cosmologies. 



1 . Introduction 



Q\ Having to accept that the expansion of the Universe is acceler- 
00 ating today and that only ~ 30% of its content is contributed 
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Generalising this, the equation of state is modified 
2 with w < -1/3. In the simplest of these models, 



by matter, we need to search for what may be driving the accel- 
erated expansion. Friedmann's equations require the dominant 
form of matter to have a pressure p < -pc 2 /3, where p is its den- 
sity and c is the speed of light. The cosmological constant has 

P = -P c? 
_^ to p = wpc 

O w is constant, but it is more natural to assume that w is a func- 
i tion of time, scale factor or redshift. One possible, admittedly 
O hypothetical form of matter with such an equation of state is a 
£j self-interacting scalar field with an interaction potential which is 
^ sufficiently larger than it s kinetic energy (e.g. [ Wetterich, 1988; 
^ iRatra and PeeblesUl988HPeebles and Ratrall2002l) 

Replacing the cosmological constant by such a hypotheti- 
cal "dark energy" has consequences for structure growth and 
, the properties of dark-matter haloes (Bartelmann et al., 2002; 
^j lWei nberg and Kamionkowski. 2003; Klvpinetal., 2003]). We 



report here on our studies of how h alo con centrations change 
in a variety of dark-energy models ( Dola g"et all I2004T) . This 
leads us to suggest a remarkably simple scaling of halo con- 
centrations with the linear growth factor in dark-energy models. 
We indicate consequences for strong lensing by galaxy clusters, 
which offer one possibility for constraining dark-energy models. 
Throughout, we use present-day matter and dark-energy density 
parameters of £2 m o = 0.3 and £2qo = 0.7, and a Hubble constant 
of h = 0.7 in units of 1 00 km s" 1 Mpc~ 1 . 

2. Dark-Energy Models 

The continuity equation requires that the density of dark energy 
changes with the scale factor a as 



£2 Q (a) = Qq exp 



J a 



[1 + w(a')]dW 



(1) 



This term replaces the usual cosmological-constant term in 
Friedmann's equation. For w = const. = — 1, the cosmological- 
constant behaviour is retained. For w = const., the scale-factor 
dependence simplifies to 



If w = -1/3, the model behaves like an open CDM model with- 
out cosmological constant because then Qq mimics the curva- 
ture term in Friedmann's equation. 

We use one model with w = -0.6 = const, for reference, and 
two m odels with time-varying w . These are the Ratra-Peebles 
model (Ratra and Peebles, 1988) in which t he scalar-field po- 
tentia l is a power law, and the SUGRA model ( Brax and Martin, 
2000) which has an additional exponential factor in the poten- 
tial. Both are normalised such that wo = -0.83 at a = 1. While 
w is almost constant for the Ratra-Peebles model in the relevant 
redshift range, it increases from -0.83 to ~ -0.4 between red- 
shifts and 2 in the SUGRA model. 

If normalised to its amplitude today, structure grows earlier in 
dark-energy compared to cosmological constant models. Since 
numerical simulations demonstrate that dark-matter haloes keep 
a memory in t heir cores of the mea n cosmic density at their for- 
mation times ( Nav arro et all fl997l) . haloes forming earlier are 
expected to have higher core densities. Thus, at fixed mass 
and redshift, haloes are expected to be more concentrated in 
dark-energy than in cosmological-constant models. This was 
expected from analytic considerations. The work reported here 
aims at testing this expectation with numerical simulations. 

3. Numerical Simulations 



Using the Gadget code ( Spr ingel et all 120011) . we ran a large- 
scale cosmological simulation of the ACDM model, identified 
massive haloes within it, identified their Lagrangian volumes at 
the initial redshift, added small-sca le power and re-sim ulated 
them at much increased resolution llTormen et all H99l . The 
particle mass is 5 x 10 9 h~ x M Q . We thus created a sample of 17 
clusters with final masses between 3 x 10 14 and 2 x 10 15 /T 1 M Q . 
The ACDM simulation was normalised to <x 8 = 0.9. 

We then re-simulated this same cluster sample in dark-energy 
cosmologies. For doing so, we shifted the initial redshift of the 
simulation to higher values such that the earlier structure growth 
in the dark-energy models was compensated. For achieving the 
same density-fluctuation normalisation today, the initial redshift 
Zini of the dark-energy simulations needs to satisfy 



D+(Zim) _ Z) +,ACDm(£™cDm) 



(3) 



Qq(<z) = Qooa 



-3(l+w) 



(2) 



D + (0) £> + ,acdm(0) 
where D + (z) is the linear growth factor as a function of redshift, 



and z™ CDM is the initial redshift of the ACDM simulation. It is 
important to also rescale the initial velocities of the simulation 
particles to the higher initial redshift. 

The normalisation of the models is an open issue. The ear- 
lier structure growth in dark-energy models increases the Inte- 
grated Sachs-Wolfe effect and thus the amplitude of large-scale 
secondary fluctuations in the CMB. A smaller fraction of the 
observed fluctuations can then be attributed to the primordial 
CMB, thus the normalisation of the power spectrum should be 
lowered. 

However, we are interested in the properties of dark-matter 
haloes whose linear scale is much smaller than that of the Sachs- 
Wolfe tail in the CMB power spectrum. How the large-scale am- 
plitude measured by on CMB translates to small scales depends 
sensitively on the large-scale slope of the density power spec- 
trum. We argue that weak gravitational lensing directly mea- 
sures the power-spectrum amplitude at the scales relevant here 
and should thus develop into the prime method for normalis- 
ing the power spectrum at small scales. For now, we study two 
sets of normalisations. One has constant <x 8 = 0.9 for all cos- 
mological models, the other has reduced erg such as to take the 
enhanced ISW effect into account. 

In total, we simulate our sample of 17 clusters in eight dif- 
ferent cosmologies: ACDM, Ratra-Peebles, SUGRA, and w = 
-0.6, the latter three with two different normalisations each, and 
an open-CDM model with £2a = for comparison. 



4. Results 

We fit the NFW density profile to all clusters and obtained con- 
centration parameters for all of them at 50 output redshifts be- 
tween z = 3 and today. For all these cluster snapshots, we 
also compute analytically e xpected concentratio n s according to 
the algorithms proposed by iNavarro et al.l(ll997l) : lBullock et af] 
( 20011): lEke et alJ ( 120011) . These algorithms implement in dif- 
ferent ways essentially the same idea: The core halo density is 
determined by the mean cosmological density at the halo forma- 
tion redshift, which is typically defined as the redshift when the 
most massive progenitor of the final halo reaches a certain small 
fraction of the final halo mass. The factor between the mean 
cosmological density and the halo core density, and the fraction 
of the final halo mass used for defining the halo formatio n red- 
shift, are tw o free parameters in t he algorithms by Navarro et alJ 
( 1997| and lBullock etalJ ( 120011) . The algorithm bv lEke et alJ 
(200 If has only one free parameter. We determine these param- 
eters by minimising the squared deviation between the analytic 
halo concentrations expected for the given halo masses and red- 
shifts, and the numerically-determined halo concentrations. 

We find that excellent agreement between the numerical 
halo concentr ations and the predictions by the algorithm of 
Bull ock et all 12001) can be achieved if haloes are assumed to 
form very early, or, in other words, if we assume that halo 
core properties are imprinted when only a very small fraction 
of the final hal o mass is already in place. The results from 
the algorithm of E ke et alJ (12001) agree very well with the nu- 
me rical concentrations w ithout any parameter adaptation, while 
the Navarro et al. ( 1997) algorithm predicts too shallow redshift 
evolution (cf. Fig.[Q. 

Our results admit a power-law fit of the form 



SUGRA, a 8 =0.9 



c(M,z) = 



M 



1 +z\l0 u h- l M e 



(4) 



where c is a constant for each cosmology. The exponent a » 
-0. 1 for all cosmologies tested, thus the mass-dependence of the 
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Fig. 1 . — Example for the agreement between analytically expected and 
numerically determined halo concentrations in the SUGRA model nor- 
malised to erg = 0.9. The solid line with the error bars shows our nu- 
merical results in eight redshift b ins betwee n redshift s 3 and 0. "NFW" 
and "ENS" stand for INavarro et alJ < 19971) and Bull ock et af] 120011) . 
respectively. 



Table 1 . — Parameters cq and a for the fit formula @ for halo 
concentrations in our eight cosmological models. 



model 


o" 8 












a 


ACDM 


0.90 


9.59 


± 


0.07 


-0.102 


± 


0.004 


Ratra-Peebles 


0.90 


10.20 


± 


0.07 


-0.094 


± 


0.005 


Ratra-Peebles 


0.82 


9.30 


± 


0.06 


-0.108 


± 


0.005 


SUGRA 


0.90 


11.15 


± 


0.09 


-0.094 


± 


0.006 


SUGRA 


0.76 


9.46 


± 


0.07 


-0.099 


± 


0.005 


w = -0.6 


0.90 


11.32 


± 


0.09 


-0.092 


± 


0.005 


w = -0.6 


0.86 


10.44 


± 


0.08 


-0.066 


± 


0.005 



halo concentrations is quite shallow for massive haloes. Tabled 
summarises the parameters we find. 

The scatter among the concentrations is large, but well- 
described by a log-normal about the mean given by 0. Quite 
independent of the cosmological mode l, the standard deviation 
of ln(c/c) is * 0.22 (cf. Fig.Q, see aisotTrn ^ll27)00l 

Interestingly, it turns out that a simple description can also be 
given for the cosmology-dependence of the parameter cq. We 
9.6, and for the other cosmological models 



find cf DM 



co 



ACDM 



£>+(Zcoll) 



D 



ACDM 



(Zcoll) 



(5) 



i.e. the concentrations scale in proportion to the linear growth 
factor, provided that the halo collapse redshift z co n is chosen high 
enough; in fact, even z co u — > oo produces good agreement with 
the numerical simulations. This reinforces the impression that 
halo properties are imprinted at very early times, when only a 
very small fraction of the final halo mass is already in place. 

The higher halo concentrations found in dark-energy cos- 
mologies are expected to have pronounced consequences for 
strong gravitational lensing by galaxy clusters. This is indicated 
by an analytic study we have carried out, and demonstrated by 
a numerical study using the same cluster sample described here 
(cf. the contribution of Meneghetti et al. to these proceedings). 
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Fig. 2. — Quite independent of cosmology, the halo concentrations fall 
on a log-normal distribution with a scatter of ss 0.22. 



5. Conclusions 

Using numerical simulations of 17 cluster-sized haloes in eight 
different cosmologies, we have tested and confirmed analytic ex- 
pectations for the dependence of halo concentrations on mod- 
els for the dark energy. Earlier structure growth in dark-energy 
models yields more concentrated haloes. If it is assumed that 
core halo properties are imprinted at very early times, when 
only a very small fraction of the final halo mass is already in 
place, the halo-concentration algorithms proposed by Bullock 
et al. and Eke et al. turn out to work remarkably well. Halo 
concentrations have a log-normal distribution with a scatter of 
« 0.22 about their mean values, quite independent of cosmol- 
ogy. For all cosmologies tested, the mean halo concentration at 
fixed mass and redshift scales in proportion to the linear growth 
factor at the halo formation time if the latter is defined to be very 
early. 
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